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We analyze a possibility of quantum criticality (gaplessness) in dimerized antiferromagnetic two- 
and three-leg spin-^ ladders. Contrary to earlier studies of these models, we examine different 
dimerization patterns in the ladder. We hnd that ladders with the columnar dimerization order 
have lower zero-temperature energies and they are always gapped. For the staggered dimerization 
order, we find the quantum critical lines, in agreement with earlier analyses. The bond mean- 
held theory we apply, demonstrates its quantitative accuracy and agrees with available numerical 
results. We conclude that unless some mechanism for locking dimerization into the energetically 
less favorable staggered configuration is provided, the dimerized ladders do not order into the phase 
where the quantum criticality occurs. 

PACS numbers: 75.10.Jm, 75.10.-b, 75.10.Pq, 73.43.Nq, 64.60.-i 



I. INTRODUCTION 

There has been a lot of interest in spin ladders for 
more than a decade by now, mainly due to their very 
intriguing critical properties. One of the most peculiar 
ones is that the existence of a gap (i.e. mass) depends on 
the number of legs. The spin excitations in a m-leg spin 
ladder are gapped if m is even, and the system is gapless 
(quantum critical) when the number of legs m is odd. For 
reviews see Ref. 00. The even-m-leg ladders provide a 
very interesting example of systems where the gap (mass) 
generation is not accompanied by a long-range order or 
apparent symmetry breaking. Such systems, known as 
spin liquids, are under enormous scrutiny. They are no- 
toriously difficult to realize in dimension more than one, 
and are believed to be relevant to the physics of high-T c 
superconductors. 3 Broadly speaking, spin ladders are in- 
teresting for studying gap (mass) generation: needless to 
say that this a very deep question in all physics. 

The particular question we address in this study is 
a possibility of "restored quantum criticality" (or gap- 
lessness) in a system (ladder) built from gapped blocks 
(dimerized chains) i' 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 Let us explain the issue 
taking a dimerized two-leg ladder as an example^ It is 
well known that a single Heisenberg spin-i chain with 
alternating spin exchange (a.k.a. the dimerized Heisen- 
berg chain) is gapped.— When two chains are coupled into 
a ladder, the system is gapped even without dimeriza- 
tion. However, quite remarkably and counterintuitively, 
the dimerized two-leg ladder can be gapless, as was con- 
jectured in Ref. [3| ■ Subsequent numerical work^i 7 -^ con- 
firmed the critical (gapless) line proposed by Martin- 
Delgado et alA Similarly, the originally conjectured crit- 
ical line in a dimerized three-leg ladder— was confirmed 
in the very recent DMRG calculations. 10 ' 11 

We need to stress one very important point: in the ear- 
lier wor 1 ' 1 1 on the intrinsically dimerized lad- 
ders, the intra-chain dimerization as well as the (stag- 
gered) dimerization ordering pattern of the whole ladder 
(cf. Figs. 11161 below) are taken for granted, i.e., as the 



model built-in assumptions. In the absence of a physi- 
cal mechanism which would lock a ladder with dimerized 
chains into a particular dimerization order, it is natural 
and physically more reasonable to consider various possi- 
ble ordering patterns on the same footing. It means that 
the possibility of quantum criticality (gaplessness) has to 
be addressed in this broader setting. A strong motivation 
comes from a recent work on the three-leg ladder cou- 
pled to phonons^ 3 - It turns out that the phonon-induced 
intrachain dimerizaton, appearing in this ladder at the 
spin-Pcierls transition, occurs into the columnar dimer- 
ization order (cf. Fig. [6] below). The latter, and not 
the staggered order, minimizes the spin-phonon Hamilto- 
nian. The columnar dimerized three-leg ladder is always 
gapped. 

So, we begin from the study of the dimerized two-leg 
ladder in Sec. [Til From analyses of the limits where the 
behavior of this ladder is known exactly, we conjecture 
that the gapped columnar dimerized phase has lower en- 
ergy. This is confirmed by the subsequent treatment 
within the bond mean- field theory^ii^ The latter ap- 
proach, when applied for the case of staggered dimer- 
ization, demonstrates a good agreement with available 
numerical results, and also yields quite simple formulas 
for analytical treatments of the problem. A similar pro- 
gram with analogous conclusions is carried out for the 
three-leg ladder in Sec. IIIII 



II. TWO-LEG LADDER 

We consider a two-leg dimerized spin ladder given by 
the Hamiltonian: 

N 

H 2 h = ^ ^ J a (n)S a (n) -S a (n+ 1) 

a=l,2 n=l 
N 

+ J ± Y,S 1 {n)-S 2 (n). (1) 
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The m-leg ladder has N rungs and m ■ N spins. We con- 
sider the situation when dimerization occurs along the 
chains (a = 1,2) only, while the coupling J± is constant 
on each rung. One needs to consider two possible dimer- 
ization patterns: alternated, when staggering occurs in 
both directions 

J a (n) = J[l + (-l) n+a S}, (2) 
and columnar, when 

J a (n) = J[l + (-l)"5]. (3) 
The two patterns are shown in Fig. [1] 
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FIG. 1: Dimerized two-leg ladder. Bold/thin/dashed lines 
represent the stronger/ weaker chain coupling J(l ±5) and 
rung coupling J±, respectively. Dimerization patterns: (a) - 
staggered; (b)- columnar. 



A. Analytical Limits 

Let us discuss the points on the (5, Ji)-plane where 
the properties of the model ([1]) are known. 

5 = 0, Jj_ = - two decoupled critical (gapless) 
Heisenberg chains. 

5 = 0, J±_ > uniform gapped ladder. 

(5=1- complete dimerization. In this limit the model 
reduces to: 

(a) the snake-like dimerized Heisenberg chain of 2N spins 
in case of the alternated staggered pattern (Fig. [3^) 

(b) a set of N/2 decoupled four-spin plaquettes in case of 
the columnar dimerization pattern (Fig. [2b) For a four- 
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FIG. 2: Completely dimerized ladder, 5 = 1. (a): Alternated 
staggering reduces the model (fl]) to a snake-like dimerized 
Heisenberg chain of 2N spins; (b): Columnar order degener- 
ates into a set of N/2 decoupled plaquettes. 

spin plaquette where its nearest-neighbor spins interact 
via the Heisenberg couplings J\ and J2, the spectrum is 
known^ The ground-state singlet has the energy 



Sa = -±l±-l yjjl + Jl-J^, (4) 

and is separated by a gap from the closest triplet state 
with with the energy E t = — Jl % J2 ■ 



As was hrst noted in Ref. 0, at the critical point 
2 J = J± = J c the snake-like chain (Fig. [5^) becomes 
just a uniform (gapless) Heisenberg model. Its ground- 
state energy is 

E HC = -2 (log 2- ~)j c N^ -0.88 J c N. (5) 

The set of N/2 plaquettes (|4]) has a lower ground-state 
energy 

Eu = -J C N. (6) 

The simple results (|5I6|) question the very existence of 
the critical line predicted in^. That line was conjec- 
tured essentially from a continuity argument and anal- 
ysis of the model in the two solvable (integrable) points 
6 = 0, J± = and 5 = 1, J± = 23 . The latter was 
assumed to be a critical (gapless) point, since the al- 
ternated staggering pattern (thus the snake-like chain 
in the limit (5=1) was assumed to be a ground state. 
The subsequent analyses addressed the issue of the crit- 
ical line (within the alternated staggering pattern) more 
quantitativelyjijiSiLSiSi while, to the best of our knowledge, 
comparisons with the energy of the columnar pattern 
were not done. As one can see from (|5l6p . the gapped 
plaquettes have lower energy than the critical (gapless) 
snake-like Heisenberg chain. Note that the decoupled 
plaquettes at 6 = 1 evolve smoothly from the columnar 
dimerization pattern at 5 < 1. 

Let us now look at the line 5 = 1, J± 7^ 2 J. It is 
convenient to parameterize couplings as 

2 J = J (l + <5 cff ), (7) 

Jl = J (l-5 e ff). (8) 

The plaquette pattern results in the ground-state energy 

En = -±J N(l + Jl + 3% s ). (9) 

The energy minimum of the snake-like configuration is 
given by the ground-state energy of the dimerized Heisen- 
berg chain. The latter is known to be equivalent to 
the massive sine-Gordon (integrable) model, perturbed 
by the marginal (logarithmic) corrections i 18 i 19 A leading 
order ansatz reads as 

1 5 4 ^ 3 

E DHC ~ -2J N(\og2- - + a -^) , (10) 
V 4 In -p- ' 

Ocff 

where from the recent numerical calculations^ a ~ 2.2 
and S ~ 110. To get a more accurate analytical expres- 
sion for -Edhc is an involved problem (for a recent anal- 
ysis and more references, seeSi). However numerically, 
the ansatz (|10[) or even the unperturbed sine-Gordon 
formula^ 

E SG « -2J a N (log 2 - 1 + 0.2728<5 e 4 / 3 ) (11) 

both work quite well. A direct inspection of Eqs. (|9I10I11|) 
(e.g., plots) clearly shows that the plaquettes again pro- 
vide the lower energy state on the line 5=1 near the 
point Ji = 2J S ef { = 0^ 
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B. Mean-Field Equations 

Away from integrable points we need to resort to ap- 
proximations. We will treat the dimerized ladders in the 
framework of the bond mean- field theory^i^ This ap- 
proach for the case of dimerized three-leg ladder is de- 
scribed in detail ini£, and the case of two legs is essen- 
tially a simplified version of the former. First, the spin 
Hamiltonian is mapped onto an interacting fermion prob- 
lem via a 2D version of the Jordan- Wigner transforma- 
tion, proposed by one usp^ Then, the phase differences 
due to hopping of the fermions around any given ele- 
mentary plaquette is approximated by 7T. The quartic 
fermionic terms c\ Ci ya c\ +1 Q Cj + i iQ are treated within the 
Hartree-Fock approximation, i.e., are decoupled using the 
single-particle (bond) parameters. The latter are defined 
as 

Q+ = ( c 2i,acl l+l a ), Q- = (c2t+l,aC 2 i+2,a)> (12) 

with a = 1,2. A single bond parameter suffices in the 
rung direction: P = (c2i,a c li a +i) ■ Fourier transforming 
along the chain direction and using the Nambu formal- 
ism, the mean-field theory results in the single-particle 
effective Hamiltonian 

H s ' c = *l.H (s/c) * fe + NC 2 , (13) 

k 

where the Hamiltonian density 7i s / c is a 4 x 4 matrix 
and the Nambu spinor = ( cf£ cfj c^f cfj ) . Here 
c" , is the Fourier transform of cL ■ To account for the 

ak la 

dimerization (doubling of the lattice spacing), the lattice 
is subdivided into two sublattices A and B. The explicit 
form of the effective Hamiltonian (fl~3f depends on the 
dimerization pattern (staggered or columnar) which is 
accounted for by an extra label (s,c) in the above equa- 
tion, and 

C 2 = J+\Q+\ 2 + J-\Q-\ 2 + J±\P\ 2 ■ (14) 

Diagonalization of TL s /° yields four energy eigenvalues 
±E s n (k),n = 1,2: 

E s 1/2 (k) = ^W + J 2 ±1 T 2Jli(Ji+ - Ji-)cos£: , (15) 

for the staggered phase. We define 

Ji± = J(1±5)(1 + 2Q±), J ±1 = J ± (1 + 2P), (16) 
W = J\ + + Jf_ - 2Ji+ Ji_ cos(2/c). (17) 

For the columnar pattern we obtain two doubly degener- 
ate energy eigenvalues ±E c (k): 

E c (k) = ~^W+j[ 1 , (18) 

The partition function of the single-particle Hamiltonian 
(|13p can be found in a closed form, and the free energy 



per spin is 

(19) 

where (3 = 1/fcsT and q is eigenvalue's degeneracy. We 
gave the above formula valid for m legs. The mean-field 
equations are derived from minimization of the free en- 
ergy (|19p with respect to the mean-field parameters Q± 
and P. These self-consistent integral equations are solved 
numerically. 



C. Analysis 

To start with, let us make the following observation: 
In the antiferomagnetic two-leg ladder ([T]) with intrin- 
sically dimerized chains we cannot think of a particu- 
lar mechanism to enforce a specific (i.e., staggered or 
columnar) dimerization pattern. So, having the dimer- 
ized chains, the ladder must admit the pattern which 
minimizes its free energy. Using Eq. (|19p . we compared 
the zero-temperature energies of two configurations in 
various parameter ranges. The characteristic results are 
given in Fig. [3] The columnar pattern corresponds the 
lower energies in all cases, so it is the thermodynamically 
stable state (phase). Now let us analyze the excitation 
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FIG. 3: Two-leg ladder: Energies of staggered and columnar 
dimerization configurations at T — 0. 

spectra of two phases. We will identify the minimal band 
gap of the eigenvalues (fT5"j) or (fT5]) with the spin gap of the 
corresponding phase. To demonstrate that this is true, it 
suffices, e.g., to calculate the uniform spin susceptibility 
from the free energy (I19|) with an external magnetic field 
added. 

We find that the columnar phase is always gapped with 
the gap A c = E c (0) given by 



A c = JsjuW + (^j) 2 p^ , (20) 
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FIG. 4: Two-leg ladder: Gaps of the staggered and columnar 
phases. 



where 



u = 1 + (Q+ 
p = 1 + 2P. 



-) + 



Q+-Q- 



(21) 
(22) 



Qualitatively, the columnar phase is similar to the uni- 
form two-leg ladder, i.e. gapped. The gap (f20|) persists 
in the limit of decoupled chains Jj_ — ► 0, as it must be, 
and disappears only together with the vanishing chain 
dimerization 8 — ► 0. If the latter limit is taken first, then 
the uniform ladder is gapped with A c = Jj_p/2. 

The staggered phase is more interesting. Its gap A s = 
El(0) reads as& 



(23) 



It vanishes on a certain critical line, even if neither 
dimeritation 8, nor rung coupling Jj_ are zero. Thus, 
the mean-field theory confirms the earlier conjecture , 
which has been so far corroborated only by numerical 
calculations^^ With the bond parameters Q±,P de- 
termined from the mean-field equations, we plot the gaps 
A c . s given by Eqs. (|20I23[) as functions of coupling ratio 
J±/J for several dimerizations in Fig. |4j Note that the 
gaps A c and A s coincide at "the decoupled chains line" 
J±/J = 0, as well as in the absence of dimerization. The 
qualitative differences between the gaps in the colum- 
nar and staggered phases are clearly seen. The former 
is always non-zero if the chains are dimerized, and in- 
creases with the growth of J±/J. The latter behaves 
non-monotonously, passing though a critical point. It is 
worth noting that the known analytical limits (discussed 
above) are recovered exactly by the mean-field equations. 
In particular, on the "decoupled" axis J± = 0, the mean- 
field parameter Q± whose physical meaning is the aver- 
age of the dimerization operator S a (n) ■ S a (n + 1), be- 
haves almost as that in the dimerized XF-chain, coin- 
ciding with the exact values of that quantity at the ends 
of the interval 8 £ [0, 1]. For the dimerized XF-chain 



FIG. 5: Two-leg ladder: Critical line where the gap of the 
staggered phase vanishes. Along with the present bond mean- 
field theory predictions, the numerical diagonalization data 
from Ref. are shown. 



the uniform term t and the dimerization susceptibility n 
at T = are given in terms of the elliptic functions and 
can be found, e.g., in Ref. In particular, 



5 = 0: Q± = ~, (25) 
5 = 1: Q+ = \, Q- =0. (26) 



For instance, at 5 = 1 (J± = 0), the ladder reduces to 
the set of decoupled dimers with the Heisenberg exchange 
2 J. Our Eqs. (I20M26|) result in the exact value of the 
gap A = 2J for this case. Similar, at the integrable 
quantum critical point 5 = 1, Ji = 2 J (uniform snake- 



like chain), we have Q+ = P, and the mean-field (|23|) 
correctly predicts a vanishing gap. 

We plot in Fig. [5] the mean- field critical line given by 



J± 2u c 



(27) 



t±8- 



V 



(24) 



along with the numerical results of Ref. 1. As one can 
see, the mean field works quite well quantitatively. It 
is surprising, but the mean field does much better than, 
e.g., the non- linear sigma model (NLcrM). The latter fails 
to predict the correct critical line.- Analysis of the topo- 
logical term of the NLcrM predicts that the line ends at 
8 = |, J± = 0. The dotted line in Fig. [5] corresponds 
to the case when in the gap equations we set the bond 
parameters to their maximal values Q± = P = \ (i-e., 
u = p = 2), so J±/J = 28. Such a naive "uniform 
isotropic limit" simplifies drastically the formulas, and 
provides a decent approximation. 
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III. THREE-LEG LADDER 

A three-leg dimerized ladder is defined by the Hamil- 
tonian: 

N 

#3L = ^2 ^ J a (n)S a (n) ■ S a {n + 1) 

a=l,2 n=l 
iV 

+ Ji^ [Sr(n) • S 2 (n) + S 2 (n) • S 3 (n)] . (28) 

n=l 

Similar to the two-leg case, we consider dimerizations 
of the whole ladder (|2I3[) which follow the staggered or 
columnar patterns, as shown in Fig. [6] 

f f 1 r- f > » * 

i i i i i i i i 

i i i i i i i i 

-t f r *- f f f r- 

| i i i i i i i 

i i i i i [ i i 

(a) (b) 

FIG. 6: Dimerized three-leg ladder. Line notations are the 
same as in Fig. [1] Dimerization patterns: (a) - staggered; 
(b)- columnar. 

Contrary to the previous case, now we don't have lim- 
its where the behavior of this ladder is known exactly. 
Although the non-dimerized limit 5 — is known to 
be gapless, and the three-leg ladder, roughly speaking, 
can be described as a "renormalized" spin-i Heisenberg 

I 



K /C (k) = ^ 



where 

Y s = (Ji+J[ + + Jx-J[S) cos(2fc) - (J 1+ J[_ + Jx-J' 1+ ), 
Y c = (Ji+J(_ + Ji-J[+) cos{2k) - (Ji+J[ + + Ji-J'J), 

and J' l± and W are given by Eqs. (|16ll7p where Q± — > 
Q'±. The mean- field equations are obtained via mini- 
mization of the free energy, given by Eq. (|19p with m — 3. 
Our analysis here follows the steps of the above two-leg 
case, so we will be brief. 

Similar to the two-leg case, the columnar order cor- 
responds to a state with lower energy of the dimerized 
three-leg ladder, thus this is a thermodynamically stable 
state. We present several plots of the zero-temperature 
energies for the two types of dimerization order in Fig. [7] 
The columnar phase of the three-leg ladder is always 
gapped when 5 ^ 0. The gap A c = £^(0) is given by 

A c = Ju8 , (32) 



chainiiii^ One can easily see from Fig. [6] that in the to- 
tally dimerized limit 5=1 the columnar dimerized ladder 
becomes a set of N/2 decoupled 6-spin plaquettes, while 
the staggered pattern evolves into a "decorated" chain. 
The latter, to the best of our knowledge, has been not 
analyzed earlier. 

We treat the problem in the bond mean-field approach 
as above. In fact, we will use directly our earlier results 
for the three-leg ladder coupled to phonons^, with mi- 
nor modifications. Along with one couple of the bond 
parameters Q±, defined by (|12[) with a = 1,3, one 
needs an extra couple of parameters Q'± for the chain 
in the middle, defined by the same equation with a = 2. 
The effective Hamiltonian density W /c (cf. Eq. (13|)) 
is now a 6 x 6 matrix, and the Nambu spinor ^\ = 
( 41 c ik 41 41 41 41 ) ■ The constant term 

Cz = J+\Q+\ 2 + J-\Q-\ 2 + \j-\Q'~\ 2 

+ l -J + \Q' + \ 2 + 2J 1 _\P\ 2 . (29) 

Diagonalization of 7i s ^° results in six energy eigenvalues 
±E s ^ c (k), j — 1, 2, 3 for each configuration: 

El(k) = E^k) = ±VW, (30) 

and 



n = 2,3, (31) 

I 

In the limit S — > 0, this mean-field approach correctly 
predicts gaplessness^, in accordance with the general re- 
sults for odd-m-leg ladders^ We should however point 
out that the gapless limit A c — > when 6 — * is re- 
covered only if we set (Q! — Q + ) — > 0. The latter 
limit, obviously correct from physical point of view [cf. 
Eq. (f24|) ]. is not automatically fulfilled in the bond mean- 
field equations. This is an artefact of the particular bond 
mean-field decoupling scheme, and is not intrinsic for any 
mean- field decoupling, as one see, e.g., from the (classi- 
cal) mean-field results of Ref. Hfl It is worth noting that 
the simple analytical approximations we derive, are free 
of these flaws. 

The staggered phase is gapped everywhere, except for a 
critical line of vanishing gap, when 6 0, J± =/= 0. Thus, 
we confirm the original conjecture^, which has been cor- 
roborated by the recent DMRG calculations of Almeida 
et aiii. The gap of the staggered phase, found from 



(-l) n J(W - W 1 ) 2 + SJl^W + W' + 2Y S / C ) + W + W' + 



4Jf 



6 




FIG. 7: Three-leg ladder: Energies of staggered and columnar 
dimerization configurations at T — 0. 




where u' is given by Eq. (l2"T]) with replacement Q± — > 
Q' ± . The columnar gap A c and the (critical) staggered 
gap A s , directly calculated from the mean- field equa- 
tions, are plotted in Fig. [H There is one subtlety related 
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FIG. 8: Three-leg ladder: Gaps of the staggered and columnar 
phases. 

to the staggered gap, absent in case of the two- leg ladder. 
Recall that we defined gap in a given phase as a mini- 
mal band gap of the energy eigenvalues of the effective 
single-particle Hamiltonian. What happens here is that 
the staggered gap 1(55]) is determined from the eigenvalue 
£| which has minimal band gap at small Jj_/J. However 
with the increase of Jj_ / J it reaches a certain value where 
"a level crossing" occurs, i.e., the band gaps of _E| and Ef 
are equal. Then, for bigger J±/J, the staggered gap is de- 
termined by the (minimal) band gap of E\ . In addition, 
since E\ = El ([50]) , the staggered and columnar gaps are 
equal, and according to ([32]). A> = A c = JuS in that 
region. In fact, however, one can see differences between 
A> and A c . The reason is that the self-consistently de- 
termined mean- field parameters entering Eq. (|21[) for u, 



are not necessarily equal for two configurations. 

To get some analytically tractable formulas for the gap, 
let us approximate u — u' . Then 



1 V2 J 1 



(34) 



Making even a more drastic approximation for the bond 
mean- field parameters by their maximal values, i.e., u — 
2 and p = 1 + l/\/2 {P = 2~§), we get 



a: 



2J\5 



1 + V2 Jj_ 



(35) 



4 J 1 

We plot in Fig. [9] the critical line obtained from the "ex- 
act" mean-field gap equation (|3"3"|) and the DMRG result 
of Ref. [U The dotted line J ± /J = 1.65695 is the pre- 
diction of the simplified gap equation (|33|) . Note that 



J/J 1 




0.2 0.4 0.6 0.8 

8 

FIG. 9: Three-leg ladder: Critical line of the vanishing stag- 
gered gap. Our results are shown along with the DMRG data 
from Ref. [ill- 
contrary to the two-leg ladder, where we have exact pre- 
dictions for the two integrable end points of the critical 
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line, for this case we don't have an independent predic- 
tion for the critical ratio Jj_/ J at 5 = We can point 
out however a very good agreement between the mean 
field and DMRG. The former via Eq. ([33]) yields 1.67 [or 
1.66 via ([35]) ], while the latter yields 1.60. 

IV. CONCLUSIONS 

We analyze the possibility of quantum criticality (gap- 
lessness) in dimerized two- and three-leg ladders. Con- 
trary to earlier studies of these models, we do not imply a 
particular ladder's dimerization pattern. In this work we 
restrict ourselves to the zero-temperature properties. Wc 
find that for a given intrinsic intrachain dimerization, the 
ladder has lower zero-temperature energy for the colum- 
nar dimerization order. This is true for both types of 
ladders (two-, three-leg) considered. For the columnar 
dimerization pattern the ladders are always gapped, i.e., 
there are no lines of quantum criticality on the (6, Jj_) 
plane. 

For the staggered dimerization order, we find that the 
ladders possess the quantum critical lines, in agreement 
with earlier analyses of this problemi^&^&^ i 10 ' 11 The 
mean-field theory we apply in this study demonstrates its 
quantitative accuracy. For the two- leg ladder, the mean- 
field critical line passes through both integrable quantum 



critical points (0, 0) and (1, 2 J) on the (5, J±) plane, and 
demonstrates a good agreement with the numerical diag- 
onalization results^ in between. For the case of three-leg 
ladder, the mean-field critical line passes through (0, 0) 
(the only integrable quantum critical point) and agrees 
well with the critical line determined through the DMRG 
calculations 

Thus, from a more practical point of view, the very 
possibility of quantum criticality in dimerized ladders 
hinges on some mechanism which would lock the dimer- 
ization order into a more energetically expensive stag- 
gered configuration. A realistic example of dimerization 
in a three-leg ladder through the spin-Peierls mechanism 
results in the columnar order and so, the gapful phase 
everywhere at non-zero dimerization >H At the moment 
we are not aware of any Hamiltonian with would provide 
such a mechanism for quantum criticality in dimerized 
ladders to occur. 
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